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We investigate the quark mass dependence of meson and baryon masses obtained from 2+1 flavor 
dynamical quark simulations performed by the PACS-CS Collaboration. With the use of SU(2) 
and SU(3) chiral perturbation theories up to NLO, we examine the chiral behavior of the pseu- 
doscalar meson masses and the decay constants in terms of the degenerate up-down quark mass 
ranging form 3 MeV to 24 MeV and two choices of the strange quark mass around the physical 
value. We discuss the convergence of the SU(2) and SU(3) chiral expansions and present the 
results for the low energy constants. We find that the SU(3) expansion is not convergent at NLO 
for the physical strange quark mass. The chiral behavior of the nucleon mass is also discussed 
based on the SU(2) heavy baryon chiral perturbation theory up to NNLO. Our results show that 
the expansion is well behaved only up to m\ « 0.2 GeV 2 . 
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1. Introduction 

The PACS-CS Collaboration has been performing Nf = 2 + 1 lattice QCD simulation with 
the nonperturbatively 0(a) -improved Wilson quarks and the Iwasaki gauge action at the lattice 
spacing of a = 0.09 fm on a (2.9 fm) 3 box[[l], ^, ||, |J. The DDHMC algorithm armored with several 
improvements allows simulations in the light quark mass region, where the resulting pseudoscalar 
meson mass is down to 156 MeV. In this report we focus on chiral analyses of our data based 
on chiral perturbation theory (ChPT). We apply the SU(2) and SU(3) ChPTs up to NLO to the 
pseudoscalar meson masses and the decay constants for the pion mass ranging from 156 MeV to 
411 MeV. The finite size effects are taken into account at one-loop level of the ChPTs. We also 
examine the quark mass dependence of the nucleon mass employing the SU(2) heavy baryon ChPT 
up to NNLO. The simulation details are given in separated reports[jl|, [5p. 

2. Simulation parameters 

We employ the 0(a) -improved Wilson quark action with a nonperturbative improvement co- 
efficient csw = 1.715 and the Iwasaki gauge action at j8 = 1.9, corresponding to the lattice spacing 
of a = 0.09 fm, on a 32 3 x 64 lattice. Simulation parameters are summarized in Table (T), where the 
pion mass and the quark masses are given at each combinations of fc u d and fc s . The quark masses are 
renormalized perturbatively in the MS scheme at the scale of I /a. The statistics of each ensemble 
is given by the Molecular Dynamics (MD) time. 

Table 1: Summary of pion mass and quark masses perturbatively renormalized in the MS scheme at the 
scale of I /a. 
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3. SU(2) and S(3) ChPT analyses for the pseudoscalar meson sector 

In principle our results should be compared with the predictions of the Wilson ChPT[Q]. How- 
ever, their expressions in terms of the AWI quark masses agree with those of the continuum ChPT 
up to NLO by a redefinition of the low energy constants (LEC)[jI|]. We then discuss our results 
based on the continuum ChPT. 

3.1 SU(3)ChPT 

In the continuum SU(3) ChPT the one-loop expressions for m n , rtiK, f%, fx contain six un- 
known LECs, Bq,/q,L4,L5,L^,L^. We determine these parameters by applying a simultaneous 
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fit to m 7C /2m a i,m^/{m u d + fn s },f n ,fK, where the finite size corrections are taken into account at 
one-loop level[§]. 

The results for the LECs are listed in Table || where the phenomenological estimates with 
experimental inputs[§], the RBC/UKQCD results^ and the MILC results[|] are also presented 
for comparison. We spot some discrepancies for the central values for the results of 14,5,6,8- Large 
errors make it difficult to draw any conclusions. We instead make a comparison in terms of the 
SU(2) LECs obtained by the conversion from the SU(3) ones. This is shown in Fig. [[]. Our results 
are l 3 = 3.47(11), l 4 = 4.21(11) and l 3 = 3.50(11), l 4 = 4.22(10) with and without the finite size 
corrections, respectively. Except for the MILC result for F3, all results are reasonably consistent. 

Table 2: Results for the LECs in comparison with the phenomenological estimates, the RBC/UKQCD 
results and the MILC results. 
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Figure 1: Comparison of the results for 1$ and I4. Black symbols denote the phenomenological estimates [|To[ 
|TT| ]. Blue ones are for 2 flavor results^, 13, Hp. Red closed(open) symbols represent the results for the 
SU(3) (SU(2)) ChPT analyses in 2+1 flavor dynamical simulations [g [§. 

Although we find that the SU(3) ChPT fit gives reasonable values for the LECs, the value of 
X 2 /dof is unacceptably large. In FigJ^ we plot the fit results for m^/m u d and f K . The strange quark 
mass dependence between the data at (k^, K s ) = (0.13754,0.13640) and (0.13754,0.13660) is not 
well described. This flaw is the main cause for the large # 2 /dof. 

In order to investigate the origin of the discrepancy, we compare the NLO contributions with 
the LO ones in Fig. [|. While the NLO is at most 10% of the LO for the pion mass, the ratio is much 
larger for the decay constants, e.g., for f K it rapidly increases from 10% at m u d = to 30% around 
^ud = 0.01. The situation is worse for The NLO contribution is about 40% of the LO one even 
at m U( j = 0, most of which stems from the loops containing the strange quark. These observations 
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mean that the physical strange quark mass is not small enough to be well controlled up to NLO in 
the SU(3) ChPT. 
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Figure 2: Fit results for (am n ) 2 / (am^f 1 ) (left) and / w (right). The black symbols represent the lattice 
results. The red and blue triangles denote the SU(3) fit results plotted at the measured quark masses. The 
cyan and orange ones are for the SU(2) case. The open and filled symbols distinguish the results at K s = 
0.13640 and 0.13660. The star symbols represent the extrapolated values at the physical point denoted by 
the vertical dotted line. 
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Figure 3: Ratio of the NLO contribution to the LO one in the ChPT expansions for the pion mass and the 
pseudoscalar meson decay constants. The strange quark mass is fixed at the physical value. 



3.2 SU(2) ChPT 

Instead of extending the SU(3) chiral expansions from NLO to NNLO, we employ the SU(2) 
ChPT up to NLO for further chiral analyses without increasing the data points. Our strange quark 
mass is close enough to the physical value to allow us an analytic expansion of the SU(2) LECs 
around the physical strange quark mass. 

For m n and f n we employ the SU(2) ChPT formulae where the low energy constants B and / 
are linearly expanded in terms of the strange quark mass: B = B^ +m s B^ and / = +m s fs 1 \ 
The K meson is treated as a matter field in the isospin 1/2 linear representation of the SU(2) chiral 
transformation, which is coupled to the pions in an SU(2) invariant way. This assumption leads to 
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Table 3: Cutoff, quark masses and pseudoscalar meson decay constants determined with m ni mK,mQ as 
physical inputs. Quark masses are renormalized in the MS scheme at the scale of 2 GeV. 





w/o FSE 


w/FSE 




w/o FSE 


w/FSE 


experiment 


a _1 [GeV] 


2.176 (31) 


2.176 (31) 


h 


132.6 (4.5) 


134.0 (4.2) 


130.7 ±0.1 ±0.36 


mjjf [MeV] 


2.509 (46) 


2.527 (47) 


f K 


159.2 (3.2) 


159.4 (3.1) 


159.8 ±1.4 ±0.44 


mf^tMeV] 


72.74 (78) 


72.72 (78) 


hlh 


1.201 (22) 


1.189 (20) 


1.223 (12) 


m s /m ud 


29.0 (4) 


28.8 (4) 











the following fit formulae for and 

2 —2 i o r rf, , ft 3 2Bm ud f2Bm ud \\ 

m K = m K + fi m m ud . f K = f j 1 +/j / m ud - - 167r 2y 2 ln ( ^ 2 ) J ' (3,1) 

where m 2 K and / are also linearly expanded in terms of the strange quark mass: fh 2 K = a m + y m m s 
and/ = /i 0) +m s /i 1) . 

We apply a simultaneous fit to m^/lm^,/^ and fx and an independent fit to m\. Evaluating 
the finite size effects based on the NLO formulae of the SU(2) ChPT for m n and f K , we obtain 

m^fi = 0.009345(27)GeV 2 , /= 124.8(51)MeV, h =3.23(21), l A = 4.10(20), (3.2) 

without taking into account finite size effects, and 

,71^ = 0.009332 (26)GeV 2 , /= 126.4(47)MeV, l 3 =3.14(23), f 4 =4.09(19), (3.3) 

including finite size effects, where is the up-down quark mass extrapolated to the physical 
point. The results for and I4 are plotted in Fig. [T] for comparison. All the lattice results reside in 
the range 3.0</ 3 <3.5 and 4.0</ 4 <4.5 except the MILC result for J 3 . 

The fit results for m 2 r /2m U( j and f n are plotted in FigQ In the SU(2) case we do not observe 
any discrepancy around am ud = 0.01. The resulting # 2 /dof are 0.43(77) and 0.33(68) with and 
without the finite size corrections, respectively. These numbers are an order of magnitude smaller 
than the SU(3) case. In Fig.^ we illustrate the relative magnitude of the NLO contribution to the 
LO one in comparison with the SU(3) case. The convergences of the SU(2) chiral expansions for 
m n and f K are clearly better than the SU(3) case. 

3.3 Results for physical quark masses and pseudoscalar decay constants 

The up-down and the strange quark masses and the lattice cutoff are determined with the 
choice of m K ,niK,mQ as physical inputs. We employ the SU(2) ChPT for the chiral analysis on m K 
and niK as discussed in the above subsection. For the Q. baryon we use a simple linear fit formula 
m = a + fim ud + ym s . 

In Table | we summarize the results for the quark masses, the lattice cutoff and the pseu- 
doscalar meson decay constants together with the experimental values. Our results for the quark 
masses are smaller than the estimates obtained by recent 2+1 flavor lattice QCD simulations^, gjp. 
We note, however, that we employ the perturbative renormalization factors at one-loop level which 
should contain an uncertainty. For the decay constants we observe a good consistency between our 
results and the experimental values within the error of 2 — 3%. Here also one-loop perturbative 
results for the renormalization factors are used. We note that work is in progress to determine the 



renormalization factors non-perturbatively using the Shcrodinger functional [15] 
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4. SU(2) ChPT analysis for nucleon mass 

For the nucleon mass we employ the SU(2) heavy baryon ChPT formula upto NNLO[p^]: 

£„2 



2 " l 7t 



+ 



327T/2 

64n 2 f 2 \m 2 J 32n 2 f 2 \m J jU 



(4.1) 



which contains six new low energy constants mo,ci ,C2,C3,gA,ei. The 0{m n ) term which is ob- 



tained by the relativistic baryon ChPT[17] gives only a small contribution so that the results are 
little affected in the following analyses. 

Since the number of our lattice data points is not sufficient for a full determination of the LECs 
by the ChPT fit, wechoose tofitmo,ci,ei wbilegA,/rc,ci,C2,C3 are fixed at the phenomenologically 
viable values. Following Ref. [18] we set = 1.267, C2 = 3.2GeV _1 and consider two possibilities 
for cy. C3 = — 3.4GeV _1 as fit-A and C3 = — 4.7GeV~ 1 as fit-B. For f % we use the pion decay 
constant extrapolated at the chiral limit with the SU(2) ChPT fit. We take two fit ranges for the 
data series with fc s = 0.13640: Range-I isforO.13781 < K ud < 0.13727 and range-II for 0.13781 < 
K"ud < 0.13700. The fit results are depicted in Fig. |4| We find that the lattice results are remarkably 
well described up to nv\ = 0.5 GeV 2 . Table |] summarizes the results for the LECs together with 
the QCDSF/UKQCD results for comparison. The nucleon sigma term is also obtained through 
Own = m 2 n {dm^ / dm 2 n ). All the results are compatible within the 2a errors. 

In Fig. ^| we separately draw the contribution of the LO, NLO, NNLO and 0[m n ) terms in the 
chiral expansion. Observe that there is a drastic cancellation between the LO and NLO contribu- 
tions which are both large in magnitude, and that the NNLO contribution monotonically increases 
as m n becomes heavier. The convergence of chiral expansion in the SU(2) heavy baryon ChPT is 
hardly controlled beyond nr % ~ 0.2 GeV 2 . 





range-I 
fit-A fit-B 


range-II 
fit-A fit-B 


mo [GeV] 
c\ [GeV- 1 ] 
ei(lGeV) [GeV" 3 ] 


0.880 (50) 0.855 (47) 
-1.00(10) -1.19(10) 
3.7(1.4) 4.2(1.4) 


0.850 (27) 0.795 (27) 
-1.08 (4) -1.34(4) 
2.9 (4) 2.4 (3) 


^ 2 /dof 


0.1(0.9) 0.0(0.5) 


0.3(0.9) 1.2(1.6) 


a Nn [MeV] 


51.4(7.6) 60.1(7.3) 


56.8 (2.7) 70.8 (2.6) 



QCDSF-UKQCD 



fit-A 



fit-B 



0.89 (6) 0.76 (6) 
-0.93 (5) -1.25 (5) 
2.8 (0.4) 1.7 (0.5) 



Table 4: Results for mo, c\,e\ together with the QCDSF/UKQCD results. 
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Figure 4: Fit results of the nucleon mass using the SU(2) heavy baryon ChPT formula (left) and each 
contribution of the LO, NLO, NNLO and 0(m%) terms in the case of fit-A with range-I (right). 
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